Observation of power-law scaling for phase transitions in linear trapped ion crystals 
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We report an experimental confirmation of the power-law 
relationship between the critical anisotropy parameter and ion 
number for the linear-to-zigzag phase transition in an ionic 
crystal. Our experiment uses laser cooled calcium ions con- 
fined in a linear radio-frequency trap. Measurements for up 
to 10 ions are in good agreement with theoretical and numeric 
predictions. Implications on an upper limit to the size of data 
registers in ion trap quantum computers are discussed. 

PACS numbers:32.80.Pj, 64.60.-i, 03.67.Lx, 52.25.Wz 

Ions confined in linear radio-frequency traps, and 
cooled by laser radiation, will condense into a crystalline 
state. Such crystals are the most rarefied form of con- 
densed matter known [Q . Besides being of inherent sci- 
entific interest for this reason, cold trapped ions have 
a growing number of applications, notably spectroscopy 
[2-pj, frequency standards and quantum computing 
[ 3,H. The existence of different kinds of phase transitions 
of these crystals has been known for some time |^|| and 
has been the subject of various theoretical and numeric 
studies [^oPJllll. Previous experimental work identified 
different crystal phases/configurations in a quadrupole 
ring trap (^] . Here we explicitly investigate the transition 
between two of these phases: the linear and the zigzag 
configurations. We report the first experimental confir- 
mation of one of the key theoretical/numeric predictions 
for the linear-to-zigzag transition, namely the existence 
of a power-law relating the critical anisotropy parameter 
to the number of ions in the crystal. Further, we discuss 
the usefulness of this power-law expression in determin- 
ing the ultimate size of a quantum logic register realizable 
using a single ion trap. 

The potential energy of a crystal of N identical ions 
of mass M and charge e confined in an effective three- 
dimensional harmonic potential is 
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where r„ = (x n ,y n , z n ) is the postion vector of the n — th 
ion, and v x , v y and v z are the trapping potential fre- 




FIG. 1. Trap electrodes include four 1-mm diameter rods 
and two endcaps. Rf trapping voltages Vo range from 60-200 
V, producing radial frequencies v r — 200 — 700 kHz. End- 
cap voltages range from 10-200 V, producing axial frequencies 
v z = 80 — 390 kHz. Outer support rods (not shown) double as 
compensation electrodes for moving ions radially within the 
trap. The insert shows a level diagram for 40 Ca + . 

quencies in the three directions. Assume the trapping 
potentials are approximately equal in the two transverse 
directions (x and y), so that v x i=s v y ~ v r , and the trap- 
ping potential in the axial (z) direction is different than 
in the other two directions. This anisotropy is character- 
ized by the parameter a = (v z /v r ) 2 . Work by Schiffcr 
jlj and Dubin [jll| predicts that the ions undergo a phase 
transition from a linear to a zigzag configuration at a 
critical anisotropy value a cr it- The predicted power-law 
scaling for a cr it versus the number of trapped ions N is: 
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where the constants found empirically by Schiffer are 
c = 2.53 and = -1.73 @. Later in this paper we pro- 
vide a simple alternative approach which produces qual- 
itative agreement with Schiffer's results but is based on 
a stability analysis of the transverse oscillatory modes, 
rather than on numerical simulations of the equilibrium 
crystalline configurations. 

We studied the zigzag phase transition using strings 
of 40 Ca + ions confined in a linear radio- frequency 
quadrupole trap. Trap electrodes (see Fig. |l] or Ref. [[hS]) 
consist of four rods and two conical endcaps. An rf drive 
at J1/27T = 6 MHz is applied to two diagonally opposite 
rods to produce a radial confining pseudo-potential |l4|] 



1 



with frequency v r proportional to the rf trapping voltage 
Vq (for example, [T^]). Axial confinement is provided by 
a static potential applied to the two endcap electrodes, 
resulting in a harmonic potential with frequency v z pro- 
portional to the square root of the endcap voltage. An- 
harmonic contributions to the axial potential energy are 
estimated to be less than 3 x 10 -7 the harmonic contri- 
bution, and are thus ignored. 

Calcium ions are introduced into the trap by intersect- 
ing an atomic beam with an electron beam. Laser cool- 
ing on the Si/2 to P1/2 transition at 397 nm (see Fig. [I]) 
cools the ions so that a crystallized string is formed along 
the trap axis. The resulting 397 nm ion fluorescence is 
detected in the horizontal plane along a direction per- 
pendicular to the trap axis by both a CCD camera and 
a photomultiplier tube. An 866 nm laser returns ions 
falling into the D 3 / 2 metastable state to the cooling cy- 
cle. The 397 nm light is produced by doubling 794 nm 
diode-laser light which is intensified by a tapered ampli- 
fier; the 866 nm beam is also produced by a diode laser. 

For each measurement of a cr it, the rf trapping poten- 
tial is lowered, keeping the endcap voltage constant, until 
the transition to a zigzag pattern is observed. The rf po- 
tential is then repeatedly raised and lowered to determine 
the reproduceability of the transition point and to rule 
out hysteresis. At the identified critical rf voltage, we 
measure the radial and axial resonant frequencies. The 
measurements are repeated for a range of endcap voltages 
and then for different numbers of ions N. 

To determine the rf trapping voltage at which the lin- 
ear crystal becomes zigzag, the voltage is lowered until 
the equilibrium position of at least one of the ions moves 
visibly (0.5-1.5 pixels « 0.3 — 1.0 (im) on the CCD cam- 
era image. Axial shifts are detected more easily than 
radial shifts; therefore, movement in the axial direction 
is used as the diagnostic to identify the phase transition. 
Moreover, observations show the spacing between the two 
innermost ions to be most sensitive to this axial reorga- 
nization. The critical rf voltage is determined to within 
±0.1 V out of 3-8 V on the synthesizer generating the rf 
trapping voltage. This 0.1 V resolution, corresponding 
to 8-16 kHz uncertainty in radial frequency, is limited by 
the voltage change necessary to move an ion visibly on 
the CCD camera rather than by the synthesizer itself. 

The radial and axial frequencies are individually mea- 
sured by applying an external drive and observing melt- 
ing of the crystal into a diffuse but still stable cloud when 
the drive frequency is resonant. Note that although the 
ions are driven at a center-of-mass frequency, coupling 
to higher order modes heats the string until it melts. 
The external drive is the output of a function genera- 
tor, attenuated and capacitively coupled through 1000 
pF to one endcap electrode. Resonant frequencies are 
measured reproduceably to within ±1-2 kHz by stepping 
repeatedly back and forth through the typical 1-4 kHz 
range in which melting is observed. Since the azimuthal 
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FIG. 2. Measured radial versus axial frequencies (points), 
at the onset of the zigzag instability, agree well with the pre- 
diction of our theoretical analysis (lines). Measurements were 
taken on seven different days, with some days dedicated to a 
particular length ion string and other days spent studying 
up to six different length strings. Theory lines pass through 
the origin and have slopes v r jv z = (a C rit)~ 1,/2 , increasing 
with ion number for N — 3 (minimum slope shown) through 
N — 10 (maximum slope shown). Error bars (see text) are 
dominated by the uncertainty in determining the critical rf 
voltage for zigzag onset while the axial frequency is held fixed. 

symmetry of the trap is not perfect, we measure two sep- 
arate radial frequencies, differing from each other by 1- 
2%. The data presented includes only the smaller radial 
frequency (corresponding to the weaker axis of the radial 
well) which is responsible for the zigzag instability onset. 

Figure |^ shows the radial versus axial frequencies at 
which the linear-to-zigzag phase transition was observed 
for string lengths up to N — 10. Figure |^ shows the av- 
erage a cr u as a function of N, and demonstrates the pre- 
dicted scaling behavior for these critical parameters. In 
both figures, the solid lines are results from our theoret- 
ical analysis (presented below) with no free parameters. 
Although the measurements of a cr u lie slightly above 
theory, the overall agreement within the 5-6% error bars 
is quite remarkable for a second order phase transition. 

In Fig. ||, vertical error bars are dominated by the 8- 
16 kHz resolution in determining the onset of the zigzag 
mode as described above, but also include (in quadra- 
ture) the 1-2 kHz measurement resolution of the ra- 
dial frequency. Horizontal errors representing the axial 
frequency measurement resolution are smaller than the 
point size and not shown, but are included when calcu- 
lating the error for a measured a cr n. In Fig. |^, each 
plotted value is a weighted average of measured a cr it s - 
Corresponding error bars are taken to be weighted aver- 
age errors rather than quadrature sums of the individual 
Ccrit errors, because they represent detection resolution 
rather than statistical errors and performing many mea- 
surements would not be expected to reduce their size. 

An unaccounted source of error arises from the thresh- 
old axial shift of 0.5-1.5 pixels required to detect a phase 
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FIG. 3. Measured values (points) of the critical anisotropy 
parameter a cr u versus ion number N agree well with a fit to 
the predictions of our theoretical analysis (line) . Each point 
is a weighted average of measured values from Fig. as is 
each error bar. 



transition. Any smaller axial shifts due to onset of the 
zigzag mode at higher rf voltages go undetected. Inves- 
tigations tracking the axial position of an ion near the 
center of a linear crystal as a function of rf trapping 
voltage were performed to estimate the size of this de- 
tection threshold error. Examples such as those shown 
in Fig. |i| reveal that once a phase transition is detected, 
the axial equilibrium moves roughly linearly with rf volt- 
age at a rate of 0.3-1.1 pixels per 0.1 V step of the drive 
synthesizer. Not surprisingly, the weaker trap settings 
- low endcap voltage and small ion number - produce 
the larger amount of movement per voltage step but also 
require the larger threshold shift. Underestimates in the 
critical trapping voltage due to this effect, therefore, end 
up of the same order as the 0.1 V resolution error already 
assigned. This detection threshold error accounts in part 
for the slightly high trend exhibited by the measured 
a cr it's. Alternatively, ignoring the detection threshold 
error, our measurement can instead be interpreted as an 
upper bound on a. crit . 

Trapping and cooling parameters were varied outside 
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FIG. 4. Axial position as a function of rf voltage for an ion 
nearest to, but not at the center of, a 4 ion string (circles) 
and a 7 ion string (triangles). Dashed lines show the steady 
state axial position for high rf voltages. Solid lines are fits to 
the points whose rf voltages lie at or below the critical value. 



their normal ranges to rule out systematic effects. Detun- 
ing the 397 nm laser as far as possible to longer or shorter 
wavelengths, while still maintaining crystallization, al- 
tered neither the measured resonant frequencies nor the 
measured critical rf voltage. Similar tests involving the 
866 nm laser wavelength showed no effects. Finally, for 
each new endcap voltage, care was taken to move the ion 
string to the radial center of the rf quadrupole potential. 
Ions not centered in this potential experience micromo- 
tion at the 6 MHz trap frequency. This driven motion, 
detected as a modulation on the ion fluorescence inten- 
sity, was used as a diagnostic to find the optimal position 
along one radial direction ju| . The diagnostic of micro- 
motion for the orthogonal radial direction was the obser- 
vation of a radial shift in ion position upon weakening or 
strengthening the pseudo-potential well, because an ion's 
position in the pseudo-potential will depend on the po- 
tential well strength unless the ion is centered Ions 
were moved radially by applying voltages to the com- 
pensation electrodes, which were typically tuned to ±1 
V; tests detuning them by 5-10 V revealed no change in 
measured resonant frequencies or in critical rf voltages. 

We now present a simple theoretical analysis of the on- 
set of zigzag instability. The equilibrium positions of the 
ions are determined by the condition that the potential 
be a minimum, i.e. W n U = (n = 1, . . . N). For the case 
of strong anisotropy (a -C 1), the ions are configured 
along a line in the z-direction. The solutions to the equi- 
librium equations in this case have been investigated by 
various authors (see, for example |l6]-|l8||); we denote the 
equilibrium position of the n — th ion by r„ = (0, 0, z n ). 

Small oscillations of the ions about their equilibrium 
positions are described by the Lagrangian [[l7||l9|] 
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where (£„, r] n ,^ n ) are the displacements of the n — th ion 
from its equilibrium position in the (x, y, z) directions, 
respectively. The coupling matrices A n m and B n ^ m are 
given by the formulae (following Ref. [[l7|) 
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where I = [e 2 / Ane M(2nv z ) 2 ] 1/3 is a length scale and 
Sn.m is the Kronecker delta. 

The eigenvectors and eigenvalues of the real, sym- 
metric, positive-definite matrix A n ^ m define the normal 
modes of oscillation of the ions along the z-direction. Be- 
cause of their importance to quantum computing, they 
have been studied in some detail The eigenvec- 

tors are defined by the formula J2m=i ^n,m^>m = fJ-pbn\ 
where fj, p > is the eigenvalue, bffl the normalized eigen- 
vector, and p (= 1,...N) the mode index (the modes 
being enumerated in order of increasing eigenvalue) . 

From the definition of the coupling matrix for the ra- 
dial oscillations, B n , m , we see that it has identical eigen- 
vectors to A n ^ m , but has different eigenvalues: 

x>^ = (^-f)e (5) 

m— 1 ^ ^ 

For these oscillations, the eigenvalues are no longer al- 
ways positive. Indeed, as a is increased, a critical value 
occurs for which one of the eigenvalues of the radial oscil- 
lations becomes zero. Beyond this point the radial oscil- 
lation is unstable and so is the linear configuration: this 
marks the onset of the zigzag mode. The critical value 
of a for which this occurs is given by the exact result 
a C rit{N) = 2/{hn — 1), where /ijv is the largest eigen- 
value of the matrix A n m for a given number of ions N. 
The value of /in in general must be determined by nu- 
merically diagonalizing the matrix A n-m for each value of 
N, but the expression for a crit (N) is exact. 

We can approximate a cr i t (N) as the power-law in Eq. |^ 
in order to compare it to our experimental measurements. 
Table I compares the constants c and (3 (see Eq. ^) de- 
rived from the measurements, from a fit to our theory 
(over N = 3 — 10), and from Schiffer's fit. We emphasize 
that no matter how well any model calculates a cr u for a 
specific N, the power-law deduced for a cr it versus N is 
still approximate and must be obtained from a fit of many 
calculated a cr jt's. In light of this, the three sets of coef- 
ficients are close, particularly for (3. Small discrepancies 
may be explained by the fact that the exact expression 
for a cr i t {N) is not a power-law, and that Schiffer's pre- 
dicted values were deduced from a fit of 10 a cr u values 
over the range N = 2 — 500 whereas we fit over just the 
experimental range N = 3 — 10. Fitting, instead, over 
all values in the range N = 2 — 100, yields the constants 
c = 2.88 ± 0.03 and f3 = -1.773 ± 0.003, which are more 
appropriate (as are Schiffer's) for large N applications of 
the power-law discussed below. 

Using the power-law expression, we illustrate how to 
estimate the maximum number of ions it is reasonable to 
confine along the trap axis. This gives one possible upper 
limit to the size of data registers for quantum computa- 
tion with on-axis ions in a linear trap. Rewriting Eq. |2| 
gives an upper limit N cr i t to the number of ions one can 



TABLE I. Scaling constants c and j3 from fits of the ex- 
perimental and theoretical data shown in Fig. ^ and from 
Schiffer's numerical results. Standard errors from our linear 
regressions are included. (The error in Schiffer's linear re- 
gression was not reported). An additional —5% experimental 
uncertainty in a cr n is included to account for the detection 
threshold error discussed earlier; it dominates the —0.2 un- 
certainty in c and the ±0.04 uncertainty in j3. 
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trap in the linear configuration: N cr i t = {v z / \fcv r ) 2 /P . 
Experimental limitations will dictate how large the ra- 
tio v z jv r can be made, but this expression should prove 
useful in optimizing future designs of ion traps. 

In summary, a simple theoretical analysis has produced 
the critical anisotropy parameter a cr it for the linear- 
to-zigzag phase transition. Experimental measurements 
have provided a confirmation of the predicted power- 
law scaling for a cr a versus TV which had until now been 
untested. From this scaling we obtain an expression for 
an upper limit to the size of linear data registers for quan- 
tum computation in these traps. 
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